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Abstract 

Pinsky (1996) JS] proved that the finite mass superdiffusion X 
corresponding to the semilinear operator Lu + (3u — au 2 exhibits local 
extinction if and only if A c < 0, where A c := A c (L+/3) is the generalized 
principal eigenvalue of L + (3 on M. d . For the case when A c > 0, it 
has been shown in Englander and Turaev (2000) [5] that in law the 
superdiffusion locally behaves like exp[£A c ] times a non-negative non- 
degenerate random variable, provided that the operator L + f3 — A c 
satisfies a certain spectral condition ('product-criticality'), and that a 
and /i = Xq are 'not too large'. 

In this article we will prove that the convergence in law used in the 
formulation in can actually be replaced by convergence in probabil- 
ity. Furthermore, instead of R d we will consider a general Euclidean 
domain D C R d . 

As far as the proof of our main theorem is concerned, the heavy 
analytic method of [H] is replaced by a different, simpler and more 
probabilistic one. We introduce a space-time weighted superprocess 
(7J-transformed superprocess) and use it in the proof along with some 
elementary probabilistic arguments. 
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Resume: 

Pinsky (1996) ^5] a prouve que le processus de super diffusion de masse 
finie X correspondant a l'operateur semilineaire Lu + (3u — au 2 possede la 
propriete d'extinction locale si, et seulement si, A c < 0, ou A c := A C (L + (3) 
est la valeure propre generalised de L + (5 dans Dans le cas ou A c > 0, 
et pour un operatour L + j3 — A c possedant une condition spectrale (de 
'criticalite-produit'), et pourvu que a et fi = Xq ne soient pas trop grands, 
Englander and Turaev (2000) [S] ont montre le processus se comporte locale- 
ment et en loi comme exp[iA c ] avec une constante multiplicative aleatoire 
non degeneree. 

Dans Particle present, nous montrons que la convergence en loi de |Hj 
peut etre renforcee par la convergence en probabilite. De plus, l'espace M. d 
est generalise a un domaine Euclidien quelconque. 

S'agissant de la preuve du theoreme principale la lourde methode an- 
alytique de [E] est remplacee par une approche probabiliste plus simple. 
Nous introduisons une renormalisation spatio-temporelle du superprocessus 
(' H -transformed superprocess 1 ) que nous utilisons dans la preuve combinee 
a des arguments probabilistes elementaires. 

1 Introduction 

1.1 Preparation 

We consider a superprocess which arises as the short life time and high den- 
sity diffusion limit of a branching particle system, which can be described 
as follows: in the n th approximation step each particle has mass 1/n and 
lives a random time which is exponential with mean 1/n. While a particle 
is alive, its motion is described by a diffusion process corresponding to the 
operator L. At the end of its life, the particle dies and is replaced by a ran- 
dom number of particles situated at the parent particle's final position. The 
distribution law of the number of descendants is spatially varying such that 
the mean number of descendants is 1 + ^^, while the variance is assumed to 
be 2a{x). All these mechanisms are independent of each other. The process 
is determined by the quadruple (L, f3,a; D), where L is a second order el- 
liptic operator corresponding to the underlying diffusion process on D. See 
Appendix A in Englander and Pinsky (1999) (Hj for a precise statement on 
the particle approximation. 

We start by presenting a formal description of the model considered in 
this article. For convenience we first recall the basic notation: let D C M. d 
be a domain and let B{D) denote the Borel sets of D. We write Mj{D) and 
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M C (D) for the class of finite measures resp. the class of finite measures with 
compact support on B{D). For \x G A4f(D), denote ||/j|| := /J.(D) and let 
C^(D) and C+(D) be the class of non-negative bounded continuous resp. 
non-negative continuous functions D — > K having compact support. Write 
C k,v (D) for the usual Holder spaces of index 77 G (0, 1] including derivatives 
of order k, and set C"?(D) := C°' V (D). 

We continue with the definition of the (L, /3, a; D)-superdiffusion, X. Let 
L be an elliptic operator on the domain D C W 1 of the form 

L:= Iv-aV + 6- V, (1) 

where aij,bi G C 1,V (D), i,j = l,...,d, for some 77 G (0,1], and the matrix 
a(x) := (aij(x)) is symmetric, and positive definite for all x G D. In 
addition, let a,/3 G C V (D), and assume that a is positive, and /3 is bounded 
from above. 

We now present our model. 

Definition 1 (Time-homogeneous superdiffusion) 

Let (X, , (J, G «M/(-D)) denote the (L, (3, a; D) -superdiffusion. That 
is, X is the unique A4f(D) -valued continuous (time-homogeneous) Markov 
process which satisfies, for any g G C^{D) , 

W<s* V {X t ,-g) =exp( M ,-u(-,t)), (2) 

where u is the minimal nonnegative solution to 

ut = Lu + f3u — au 2 on D x (0, 00), 

As usual, {v,f) denotes the integral J D v(dx) f(x). 

(See Dynkin (1991, 2002) 0, |lj or Dawson (1993) for the definition 
of superprocesses in general; see Englander and Pinsky (1999) 6 for more 
on the definition in the particular setting above.) 

Remark 1 (Time-inhomogeneous superdiffusion) The model under con- 
sideration is a time-homogeneous process. However, it is important to point 
out that for the formulation of the main theorem and the proof, the in- 
troduction of certain time-inhomogeneous superdiffusions is required. The 
previous definition will therefore be generalized for time-inhomogeneous su- 
perdiffusions in Appendix B. (See Definitional) 
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Let 



A c := A C (L + 15) 

:= inf{A el : 3u > satisfying (L + (3 - X)u = in D} ' ' 

denote the generalized principal eigenvalue for L + /3 on D. Let £ L be the 
diffusion process on D corresponding to L, and denote by ¥ x the law of £ L 
starting at x G D. Then from a probabilistic point of view, the generalized 
principal eigenvalue can be equivalently expressed as 



A c = sup lim - log M x 

{A: ACCD, dA is C^}*^ 00 * 



ex V [[ p(£)ds];T A >t\, (5) 
Jo J 



for any x E D, where t a = wf{t > : £ L (t) A}, and the C^^-boundary 
is defined with the help of C 2,r, -maps in the usual way. (See Section 4.4 in 
Pinsky (1995) ^1] on the subject). Hence, since (3 is bounded from above, 
A c < oo; and it is known from standard theory that for any A > A c , there 
exists a function < / G C 2 ' V (D) such that (L + 0)f = Xf on D. (See 
Section 4.3 in Pinsky (1995) HU) 

Pinsky proved that X exhibits local extinction (i.e., the support of X 
leaves any given bounded set, P^-a.s. for each /i € M. c ) if and only if 
A c < 0. (See Theorem 6 and Remark 1 in Pinsky (1996) |15j.) 

From now on we are interested in the situation where X does not exhibit 
local extinction. We therefore assume that A c > 0. We get a first rough 
impression about the local growth rate by the following statement taken 
from Theorem 7(b) in Pinsky (1996) |15j : 

Lemma (Local behavior in expectation) 

For /i £ M. C {D), and g 6 C^(D), satisfying ^ and g j£ 0, 

limsnp eM-ptm(Xt, 9 }] = \ ° (6) 
t^oo y oo ij p < A c 

(See also Appendix^) 

We are therefore going to concentrate on scaling with the exponent p = 
X c . In addition to the concept of the generalized principal eigenvalue we will 
then need some further ones from the so-called criticality theory of second 
order elliptic operators. In particular, we will use the concepts of critical and 
product- critical (or product-L 1 critical) operators. Recall that the operator 
L + (3 — A c is called critical if there exists a positive function / satisfying that 
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(L+f3 — A c )/ = but there is no (minimal positive) Green's function for the 
operator L + (3 — A c . In this case / is unique up to constant multiples and is 
called the ground state. The operator L + [3 — A c is called product-critical if 
it is critical with ground state < <j) c , and 4> c and 4> c (i.e. the ground state 
for the formal adjoint of L + (3 — A c ) satisfy (dx, 4> c <\>^) < oo. In this case we 
normalize them by (dx,(f) c (j) c ) = 1. 

If L + (3 — A c possesses a Green's function, then it is called subcritical. 

For the reader, it will be handy to have Appendix 2 of [5] at hand, where 
a review on criticality theory is given. For a complete presentation of the 
theory, the reader is referred to Chapter 4 in |14j . 

1.2 Motivation 

When D = M. d and L + — A c is product-critical, it is known from Theorem 
1 in |H] that if |[ a<p c || 0O < oo and the initial state [i is such that {fi, (ft c ) < oo, 
then the following holds in the vague topology: 

lim exp[-A c i]X t (dx) = 4> c dx in law, (7) 

where the limiting non-negative non-degenerate random variable A^^ was 
identified with the help of a certain invariant curve. 

It is important to point out that even though product-criticality is equiv- 
alent to the ergodicity of an auxiliary diffusion process (see next section), 
the original motion process corresponding to L does not have to be ergodic. 
In fact it can even be transient - see Example 23 in |Sj. 

(On the other hand, it follows from the discussion in Appendix A that 
when L + (3 — A c on D C M. d is not product-critical, then for /i G A4 C , and 
g G C+(D), 

lim exp{-\ c t](X t ,g) = in L 1 .) (8) 

t^oo 

There are two disadvantages of the method used in [Sj. First, the as- 
sumption that D = M. d is essential. Second, the proof does not yield any 
probabilistic insight. 

In this paper our goal is to improve the statement in and to provide 
a proof that is probabilistic in nature. We will show that the l Law of Large 
Numbers'' holds, that is, that one can replace the convergence in law by 
convergence in probability 1 . Furthermore we will drop the assumption that 

1 Since the limit is not constant, therefore, unlike in classical probability theory, one 
has to distinguish between convergence in law (WLLN) and convergence in probability 
(LLN). 
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D = M. d . In the proof we will replace the analytic reasoning given in jS] 
(which relies on dynamical systems) by a more probabilistic one using space- 
time weighted superprocesses (//-transforms). 

We suspect that in fact the Strong Law of Large Numbers holds, that 
is that the convergence in probability can be replaced by almost sure con- 
vergence. However we could not upgrade the proof of this paper to give the 
Strong Law. 

In the recent paper ^U] the authors study a supercritical superprocess 
taking values in the space of finite measures on [0, 1] , whose underlying 
motion is the Wright-Fisher diffusion corresponding to the operator 

1 , N d 2 

and whose branching mechanism is ju{l—u) with 7 > (that is, a = j3 = 7). 
They establish a dichotomy in the long-time behavior of this superprocess. 
For 7 < 1, the mass in the interior (0, 1) dies out after a finite random 
time, while for 7 > 1, the mass in (0, 1) grows exponentially with rate 7 — 1 
(as t — > 00 and with positive probability) and is approximately uniformly 
distributed over (0, 1). 

This result is in line with that of [S] if one considers the restriction of 
the superprocess on the (open) domain (0, 1). In fact it is easy to show that 
A c := 7 — 1 is the principal eigenvalue of the linearized elliptic operator L+7. 
Here is a possible argument: the operator L + 7 — A c = L + 1 can be h- 
transformed {h = v, where v is an explicitly given function in the paper) into 
a diffusion operator, which - according to their Lemma 20 - corresponds to a 
(positive) recurrent diffusion. Consequently, this /i-transformed operator is 
critical, and thus its principal eigenvalue is zero. By /j-transform invariance, 
the same is then true for the original operator L+7— A c . (See again Chapter 4 
in |14j). Furthermore, the product-criticality and boundedness assumptions 
are automatically satisfied by the boundedness of D = (0, 1). 

Finally the fact that the limiting measure is the Lebesgue measure, is also 
in line with jSj. Indeed, according to [jjj, the limiting density is a harmonic 
function with respect to the adjoint of L + 7 — A c = L + 1, that is with 
respect to L + 1, where L is the adjoint of L. An easy computation reveals 
that 

£ + l = ix(l-x ); ^ + (l-2x)A. 

Since the adjoint of a critical operator is also critical, and since positive har- 
monic functions for a critical operator are unique up to constant multiples, 
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the limiting density must be a properly normalized constant on the unit 
interval, that is, the limiting density is one. 

However, as the authors point out referring to [SJ, ' their methods use in 
an essential way the fact that their underlying space is M. d ( and not an open 
subset ofM. d , like (0, 1)), and therefore their results are not applicable to our 
situation. ' 

In the present article, as already mentioned, we manage to overcome this 
difficulty, so the result of ^Uj will fit our main result (in |1L)| the result is 
somewhat stronger as they prove convergence in L 2 ). 

2 Main Result 

Recall from (JIJ the definition of the principal eigenvalue A c of L + (3 on D 
and the corresponding ground state 4> c , and that throughout the paper we 
assume that A c > 0. Also, {<S s } s >o will denote the semigroup ('expectation 
semigroup') corresponding to the operator L + (3 on D. So far we have 
recalled (J2J) • In order to replace in ((7J) the convergence in law by convergence 
in probability, we will assume the same conditions as in jSj, except that we 
work with a generic domain D. 

Assumption 1 

In addition to the assumption that A c > 0, also assume that L + (3 — \ c 
is product- critical, that acj) c is bounded and that X starts in a state n 
with (fi, (fi c ) < co. 

Before reading the remainder of this section, it is recommended that the 
reader consults Appendix B regarding the definition of time-inhomogeneous 
superdiffusions as well as the space-time ii-transform (weighted superdiffu- 
sion) introduced there. 

Let A be a (L, /3, a; D)-superdiffusion with Xq = fi. Let H(x,t) := 
exp(— \ c t)(j) c (x), x £ D, t > 0. It turns out (see LemmaOJin Appendix B) 
that the (time-inhomogeneous) process X H defined by 

/ dA^ \ 

Xf := H(-, t)X t (that is, = H(,t)\ , t > (9) 

is an (L + a^^- ■ V, 0, a<p c e~ Xct ; D)-superdiffusion. In the sequel E and Var 
will denote expectation and variance with respect to the law of X H . 
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Lemma 1 (Bounded variance) 

; l>00 

lim Vax (||Xf ||) = / dse~ 2AcS (jM,S a [a<^]) < oo, (10) 

and IjX^H is a uniformly integrable P 9 ^ -martingale. 
Proof: Let X denote the total mass process, i.e., 

X H :=\\X H \\. (11) 

Abbreviate 

L fr :=L + a YA. v (12) 

(Pc 

and note that in fact 

Lf{u) = 4>-\L + (3- \ c )(<f) c u) = H-\L + (3 + d t ){Hu). 

(Here dt denotes differentiation with respect to time.) Let S^ c denote the 
/i-transformed semigroup with h = (f) c , that is Sf c (-) = (0 c ) _1 iS s ((/> c -). 

Define := e~ XcS Sf c ; then the semigroup {>S^} s >o corresponds to the 
operator Lq c that has no zeroth order part. In particular then 

Sf 1 < 1. (13) 

Finally, the product-criticality assumption on L + (3 — A c guarantees that the 
diffusion process corresponding to Lq c on D is positive recurrent (ergodic) 
(see Section 4.4. in Pinsky |14j). (Since ergodicity implies conservativeness, 
thus in fact Sfl = 1; nonetheless, for us it will be enough to know Qlcij) .') 

By Lemma O^b) of Appendix B along with Theorem A2 in jH], we have 
that for all / G C 2 onst (D) := {/ G C 2 (D) : CC D such that / = 

const on D \ £1}, 

d(Xf ,/) = (X?,4°f)dt + dM t (f), (14) 

where {Mt(f)}t>o is a square-integrable P'^-martingale, and its quadratic 
variation (i.e. the increasing process in the Doob-Meyer decomposition) 
(M(/)) is given by 

(M(f)) t = f dse~^ s {X^,a^ c f 2 ), t > 0. (15) 
Jo 
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(One can take the function class C^ onst (D) instead of just C^(D), because 
the diffusion process corresponding to Lg c on D is conservative, that is, it 
never leaves the domain D with probability one.) 

Applying ((Tl)) to the function / = 1, it follows that X is a P^- 
martingale. Furthermore, by (|15j). 

E^[(xf,l) 2 ] =0Mc> 2 + /'dse"^ <0 c /i,S s H (16) 

JO 

That is 

Va^(||Xf||)= [ t dse~ x ^(4>^,Si I [a ( j )c ]}= f'dse' 2 ^ {^S s [a<i^). 
Jo Jo 

(17) 

Letting t — > oo we obtain the first statement of the lemma. 

Replacing t by oo in the first of the integrals in (|17[) , we have from (|13|) 
and from our assumptions that 

Var (||Xf ||) < / d S e~ AcS (<£ c //,<Sf [c^ c ]} < A^ 1 || a<£ c |U <p,0 c ) < oo. 

JO 

Hence, by (fTo]) . 

sup t > E^ [(Xf,l) 2 ] <oo, 



and consequently X is uniformly integrable. This completes the proof of 
the second statement of the lemma. □ 

Remark 2 Our proof of LLN will indeed use the condition that a<p c is 
bounded, however it is quite possible that this condition is not necessary and 
that assuming the finiteness of the integral in (|1()[) (along with (/x, <p c ) < oo) 
would suffice. 

An immediate consequence of uniform integrability is that E^fX^J = 
(fi,4> c ), which is finite by assumption, and positive for /j / 0. This yields 
that P^fX^ = 0] < 1 for fi ^ 0. We record this in a lemma 

Lemma 2 (Limit of the total mass) 

The martingale X H has a P^-a 
is positive with positive probability 



The martingale X has a P^-a.s. limit I M := lirn^oo X t which 
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2.1 Heuristics for the Law of Large Numbers 

Before stating the Law of Large Numbers for the class of superdiffusions 
under consideration, in this subsection we give some heuristic computations. 
These will justify why we call our main result 'the Law of Large Numbers'. 

Working with the i^-transformed superprocess and at the same time, 
having the particle approximation in mind, consider particles with underly- 
ing motion Y corresponding to the elliptic operator Lq c (the probabilities 
for Y will be denoted by {¥ x , and with critical binary branching 

at rate exp[— X c t]a(x) at position x £ D and time t > 0. Furthermore let 
the system be started with initial discrete measure being "close" to v ^ 0. 

Let denote the collection of particles alive at time t in the n th ap- 
proximation step. Finally, the event of 'survival' is {[2^1 > Vi > 0}. 

The Law of Large Numbers would mean that if ^ / £ C^~(D), then as 
t — > oo, (and without further specifying what "~" means), 

M J timS > al 0,1 <w |>0V(>0 >- < 18 > 

Now recalling that in the n th approximating step the individual particle 
mass is scaled down by n and recalling also Lemma 121 one has that (for n 
large), |2~"| nI M as t — > oo. Putting this together with ((TH)) . one gets 
formally, that for large n, 

n * 

"F,E*[/(y t )]> 



v \\ /.. MfJ /vMi ~ X oo> ast^oo. (19) 



Note that in fact 

(u,W[f(Y t )]) = E,"(Xf, f) = e-^B^XtJ^} 

(y = c /u). (The first equality can be shown for instance by taking first 
the particular case v = S x , x £ D, and using that the two expectations 
satisfy the same parabolic problem; then integrating with respect to v(dx).) 
Furthermore, passing to the limit (as n — > oo) formally, the numerator of 
the fraction on the left hand side of ((191) becomes 

(X t H J)=e- x c t {X t Jcj )c ). 

Hence, for the new test function ^ / := f(j> c £ C£(D), 



(X t ,f) _x 
w(x t J) " 



H 
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2.2 Main theorem 

Making the intuition of the previous subsection precise, we now state our 
main result: 

Theorem 1 (Law of Large Numbers) 

Let f G C+(D). Iff^O and \\fi\\ ^ 0, then 

lim = in P^-probability. (20) 

Comparing our theorem with l|7|). we can now identify the limiting dis- 
tribution: N 11 = in law. 

Remark 3 

One has to be a bit careful though when making heuristic inferences 
using the particle picture as in the previous subsection. 

Obviously, the discrete system in the n th approximation step is so 
that liminf^oo > 1 under survival. That is, liminf|_ >00 \ \I-t \ > 

under survival. Recall that, heuristically, (for n large), ^;\1t \ ra -^oo 
as t — > oo. 



On the other hand, in the recent paper [5] an example of a superprocess 
is given that satisfies the conditions of our previous theorem and for 
which 

P»(xl = | 5) > 0, 
where S is the event of survival, S := {\\X t \\ > 0, Vt > 0}. O 

Conjecture 1 

We conjecture that convergence in probability in (|2T)|) can be replaced 
by almost sure convergence. 

We close this section with a remark concerning an old result. 
Remark 4 

A simple case of a super diffusion is when D = M. d , d > 1, L = iA, 
with a,f3 positive constants (supercritical super-Brownian motion). 
Here A c = (3 and 

iA + /? -A c = lA. 
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Since C = (j) c = 1, d > 1, the operator 2 A is either critical but not 
product-critical {d < 2), or subcritical (d > 3). Therefore this case 
is not included in our setup. On the other hand, the correspond- 
ing (Strong) Law of large Numbers is well known - at least for the 
discrete particle systems. Using techniques from Fourier transform 
theory, Watanabe ( ^7j) proved SLLN for branching-Brownian motion 
in M. d and in certain subdomains of it. It is not clear however if his 
method can be generalized for more general branching diffusions. o 



3 Proof of the result 



The proof is based on two observations. The first one is that the problem 
can be formulated in terms of X H , that is, one can reduce the problem to 
the investigation of a critical superdiffusion with ergodic motion component 
and exponentially decaying branching rate (again, recall that X H is an (L + 
a ^jt ' V ' °> «0ce~ Ac '; ^-superdiffusion). 

The second one is that by considering some large time t + T (where 
both t and T are large), the changes of X are negligible after time t, while 
the remaining time T is still long enough to distribute the produced mass 
according to the ergodic flow given by the if-transformed migration. 

To simplify notation, we will write W := X H (and, accordingly, Woo '■= 



-H 



X^). Denote v := Wq. By assumption, 
show that for all e > 0, 



(/J,,4> c ) < oo- We need to 



lim 



(Xt,f) 



W{X t ,f) 



> e 



0. 



(21) 



Recall that P denotes the probabilities with respect to the law of W. De- 
noting f =: g G C+(-D), we rewrite (|2*T|) as 



lim V v 



(W t ,g) 



W r . 



> e 



0. 



^(W u g) IM, 

It is easy to check that the limiting invariant density for the diffusion cor- 
responding to Lq c is 0c<^c (recall the normalization f D dx 4> c (x)(j) c (x) = 1). 
Since Z := WW is just the deterministic Lq c -Aow starting from v, therefore 

lim E v (W t ,g) = I M| • (<f> e $ c ,g), 

t— too 

and consequently our statement is tantamount to saying that for all e > 0, 



lim P l 

t— >oo 



(W t ,g)-(<t>c<l>c,g)W c 



> e 



0. 
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Let T > and let Z\y t denote the deterministic Lq-Aow starting from the 
(random) measure Wt- Then 

P u (\(W t+T ,g) - (</> c ^ c ,g)W 00 \ > e) < S x (t) + S 2 (t,T) + S 3 (t,T), (22) 

where 

Si(t) := P u (\(^c, 5)^oo - {<t>ch,9)\\Wt\\\ > e/3 
S 2 (i,T) := P"(|<^c,0>||W t || - (Z Wt (T), 5 >| > e/3 
S 3 (t,T) := P^(|(Z m (r), 5 ) - (W t+T ,g)\ > e/3). 
Take lim sup t _ >00 lim supj.^^ on both sides of Q22|). We have 



lim sup 

i— >oo 



where 



Now, since ||Wt 



(W t , 5 )-(0 c C)5 )W c 



I := lim sup Si (i) 
1/ := lim sup lim sup 52 (t,T) 
III := lim sup lim sup S3 (t,T). 



Wno as t — * 00 a.s. 



> e <I + II + III, 



(23) 
(24) 
(25) 



/ = lim Si(t) = 0. 



Also II = 0, because for all fixed t > 0, limr-»oo ^(t, T) = 0. (Indeed, 
for a/Zw G O, \\m T ^ OQ {Z Wt ^- ) (T), g) = (4> c (f> c , g)\\W t (uj)\\.) Therefore, if we 
show that 

III = 0, (26) 

then we are done. 

In order to do this, use at time t that W is a time-inhomogeneous 
Markov-process, and then apply Chebysev's inequality: 



5 3 (t,T) = &P W * (\(Z Wt (T),g) - (W t+T ,g)\ > e/3) 



<9e" 2 E^2 (W T ,g), 



(27) 



where P is the law of the (L + a^f- ■ V,0,a(f) c e Ac (*+ S ); Z?) -super diffusion 
(here t is fixed and s > is time) and a 2 denotes variance. 
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Let us now recall how the formulae for the first two moments of (Wt, <?) 
are obtained: by writing ug(t,x) for the solution of the semilinear parabolic 
evolution equation (corresponding to the superprocess) with initial value 9g, 
one differentiates (repeatedly) with respect to 8 and sets 6 = 0. 

For time-homogeneous processes with constant branching rate this is 
written down in detail in j^j, p. 37. Since the derivation of these 'moment 
formulae' only uses differentiation with respect to 9 (and not t or x), there- 
fore the proof goes through for the more general setting where coefficients 
are space-time-dependent. 

In our case, from these moment formulae and from l|27|l. one obtains 
(recall also ifTSjl ) that for all T > 0, 



(28) 



S 3 (t,T) <\-W f T &s2e- x ^{W t) SX [a<p c (S*_ s g) 2 ]) 
e Jo 

C-E u \\W t \\ 
with 

C = C(|| fl ||,||a0 c ||) := 18||a<A c || • || 5 || 2 - 

(Note that we have an extra factor 2 relative to 0| — indeed in the non- 
linear term in the semilinear parabolic evolution equation is ^u 2 .) Recall 
that ||W|| is a P^-martingale with mean and continue (|28[) with 

C- \\v\\ 



e 2 • X c e Xct ' 
Since this holds for all T > 0, thus 

limsupS 3 (i,T) < 5" 

Letting t — ► oo, one obtains (|2fjjl. completing our proof. □ 

A The behavior of the process in expectation 

For the cases when L + (3 — A c is subcritical, or critical but (dx, (f> c (t>c) = °o> 
Theorem 7(b) (ii) in ^5] states that for g £ C+, 

lim e- Xct W[(X t ,g)} = 0. (29) 
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Note, however, that in the proof there is a glitch: the proof simply refers to 
Theorem 4.9.9 in ,14 which deals with the product-critical case only, and is 
therefore not applicable for the cases mentioned. 

Nevertheless, for the subcritical case, and for fi € Ai c (D), the statement 
can be verified by a very simple argument as follows (cf. Theorem 4.9.1. in 
|14j). First, note that by the first moment formula, 

e- X ^[(X t ,g)] = { f ,, [%g](x)), 

where {%}t>o = {e~ Xct <St}t>o denotes the semigroup corresponding to the 
operator L + {3 — A c on D. 
Make an /i-transform now: 

[T t g}{x) = h(x)[T t h (g)]{x), 

where g := gh^ 1 , and pick an h > satisfying (L + j3 — X c )h = 0, to reduce 
the problem to the proof of 

]im(h^E x [g(tf)])=0, (30) 

t— »oo 

where £ h denotes the diffusion corresponding to T h , that is, to the operator 
Lq (defined analogously to Lq c in (|T2)) with (j) c replaced by h) and E denotes 
the corresponding expectation. (Of course, hfi € Ad c (D)). 

Furthermore, it is enough to show the statement with /j, := S x , x G 
D, because once we know that, the general statement follows by bounded 
convergence: E 31 ^^)] < \\g\\ for all x £ D and t > 0. 

In |14j . Chapter 4 it is shown that subcriticality is invariant under h- 
transforms, and that the transience of a diffusion is equivalent to the sub- 
criticality of the corresponding elliptic operator. Therefore, in our case, it 
follows that £ h is transient. Since g is compactly supported, by an obvious 
comparison argument, it is enough to show (J30J) with g replaced by the in- 
dicator 1b, where B is a ball. By transience 1b (£,t) — > as t — > oo a.s., and 
the statement follows from this and bounded convergence. 

Similarly, the critical but non-product critical case can be reduced to the 
analogous (but much subtler) problem of showing (|3U|) for a null-recurrent 
£ h . (Cf. the well known analogous limit theorem for countable state space 
Markov chains - see e.g. Proposition 5.3 and Corollary 6.39 in 11 .) In 
fact, as mentioned in the notes at the end of Chapter 4 in ^3], this result 
is known in the case when L is symmetric with respect to some reference 
measure pdx (see or ^Hl)- (Recall that L is symmetric if and only if 
b = aVQ for some Q € C 2,ri (D), 7] S (0, 1], and in this case L possesses a 
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self-adjoint extension due to the Friedrichs extension theorem - see Chapter 
4 in |14|.) Recently Pinchover completed the result by proving it for the 
general (non-selfadjoint) setting (see |13j). 

Consequently, the p > X c (over-scaling) part of © is immediate. One 
does not need however the above deep result for the p > A c part. Here is 
a simple alternative proof: using an h-transform with an h > satisfying 
(L + (3 — X c )h = 0, the statement is equivalent to 

lim e^-rtE'Wtf)] = 0, (31) 

t— »oo 

for each x € D, which is true in virtue of the boundedness of g. 

The p < A c (under-scaling) part of © is harder, and we are only able to 
provide the rigorous proof of the somewhat weaker assertion: 

lim sup exp[-pt]&>[(X t ,g)] = do. (32) 

t— >oo 

To this end, denote by p L+l3 (t; x, •), x 6 D, the kernel corresponding 
to the operator L + [5 and note that since g is compactly supported, by an 
obvious comparison argument, it is enough to prove that 

limsup-log / p,(dx)p L+l3 (t;x,B) = A c , (33) 

for each x € D, and Borel set B CC D. Clearly, we may assume that 
\\p\\ = 1. To verify (J33|). make again an /i-transform with an h > satisfying 
(L + — X c )h = on D. Then L + (3 transforms into Lq + A c . Moreover, 
since the generalized principal eigenvalue is invariant under /i-transforms 
and since 

A c (Lq + X c ) = X C (L^) + A c , it follows that A C (L#) = 0. Let 
PQ(t,x, •) denote the transition measures corresponding to Lq. Fix x & D 
and i? CC -D. Since h(-)h~ l (x)p L+ ^(t,x,-) = e Xct pQ(t,x, ■) (see Theorem 
4.1.1. in Jl]), and since /i is bounded between two positive constants on B, 
we have 

limsup-log J p(dx)pQ(t; x, B) = lim sup - log J p(dx)p L+ ^(t;x,B) — X c . 

(34) 

Since \\p\\ = 1 and pQ(t;x,B) < 1, the left hand side of (|3"1|) is non-positive, 
giving immediately 

X + {x,B) := limsup-log / p{dx) p L+f3 (t; x, B) < X c . (35) 
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Suppose now that X+(x, B) < A c and pick 

ce(-\ c ,-X+(x,B)). (36) 
Then by (|34[) and (J36() . along with Fubini's theorem, one obtains 

f' /*oo 

/ fx(dx) / e c V L+/3 (^,2;,5)di < oo. (37) 
Jb Jo 

Hence, for almost every x £ B, 

/•oo 

G {L+(3+c) ^ ^ ._ / e ct p L+l3^ XiB ^ dt< 00 . ( 38 ) 

Jo 

It follows from general theory then that e ct p L+ @(t, x f , B') is in fact integrable 
for all x' € D and B' CC D, that is, that the operator L + (3 + c possesses 
a (minimal positive) Green's function on D; however this contradicts to 
the well known fact that an operator with positive generalized principal 
eigenvalue does not possess a Green's function. (In our case A C (L + (3 + c) = 
A c + c> 0.) 



B The .^-transform of superdiffusions 

This section treats a generalization of the /i-transform for superdiffusions 
introduced in [B]. The /i-transform was used in the proofs in jSJ. The method 
used in the present paper however requires the spatial function < h = h(x) 
to be replaced by a space-time function < H = H(x,t). (The reader 
should not confuse with the space-time harmonic transformation yielding a 
Girsanov-type change of measure - see e.g. [T2].) 

We start with the more general definition of a time-inhomogeneous su- 
perdiffusion. Let L be an elliptic operator on D x R + of the form 

L:=Iv-aV + 6-V (39) 

where the functions a,ij,bi : D x M + — ► R, i,j = l,...,d are C 1,V (D) (for 
some rj £ (0, 1]) in the space, and continuously differentiable in the time co- 
ordinate. Moreover assume that the symmetric matrix d(x,t) := (ai,j(x,t)) 
is positive definite for all x £ D and t € M + . 

In addition, let a,/3 : D X M + — > M, be C^^D) in the space, and con- 
tinuously differentiable in the time coordinate. Finally assume that a is 
positive, and (3 is bounded from above. 
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Definition 2 (Time-inhomogeneous (L, (3, a; D)-superdiffusion) 

(i) The (L, $, a; Z))-superdiffusion is a measure- valued (inhomogeneous) 
Markov process, (X, P^ r ; p G M f (D), r > 0), that is, a family {P^' r } 
of probability measures where P^' r is a probability on C([r, oo)) and 
the family is indexed by Mf(D)) x [0, oo), such that the following 
holds: for each g G C^(D) and /i G Mf(D), 

E^[exp-(X t ,g)}=exp-(fi,u(;r;t,g)), (40) 

where u = u(-, -;t, g) is a particular non-negative solution to the back- 
ward equation 

— d r u = Lu + j3u — au 2 in D x (0, t), 
limu(-,r;i,5) = g(-). ^ 

(ii) To determine the solution u uniquely, use the equivalent forward 
equation along with the minimality of the solution: fix t > and 
introduce the 'time-reversed' operator L on D x (0, t) by 

L := iv-aV + 6- V, (42) 

where, for r G [0, t], 

a(-, r) := a(-, t — r) and &(•, r) := b(-, t — r); 

furthermore let 

/?(-, r) := /?(•, t — r) and a(-, r) := d(-, t — r). 

Consider now v, the minimal non- negative solution to the forward 
equation 

d r v = Lv + (3v - av 2 in D x (0, t), 
]imv(-,r;t,g) = g(-). ^ 

rj.0 

Then 

u(-,r;t,g) = v(-,t- r;t,g). 
(See also the remark following this definition.) 
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Remark 5 (Minimal non-negative solutions for forward equations) 

In the time-homogeneous case, minimal non-negative solutions for for- 
ward equations have been constructed in Appendix A in jHJ , and in Sec- 
tion 2 in [7] — the construction uses the approximation of D by com- 
pactly embedded subdomains with Dirichlet condition on their bound- 
aries. The construction goes through for the time-inhomogeneous set- 
ting. 

In [Hj and [7j the time interval is [0, oo) rather than [0, t]. However that 
does not make any difference — in fact the solution on the infinite time 
interval can be defined by first working on finite time intervals and then 
showing that they can be taken arbitrarily large without having the 
solution blown up. 

As we will see in Lemma|21(b), one way of defining a time-inhomogeneous 
superdiffusion is to start with a time-homogeneous one, and then to apply 
an 1 H -transform' . In general, the iif -transform of a time-inhomogeneous 
superdiffusion is defined as follows. Let < H G C 2 ' V (D) x C 1 ' V (R + ) and 
let X be a (L,(3, a; D) -superdiffusion. We define a new process X H by 

/ dX H \ 

X? := H(;t)X t (that is, = H(-,t)\ , t > 0. (44) 

This way one obtains a new superdiffusion, which, in general, is not finite 
measure-valued but only a-finite measure- valued. That is, if A4(D) denotes 
the family of all (finite or infinite) measures on D, then 

X? G M { S(D) := {v G M{D) \ H(-,ty l u G Mf(D)} 
(c.f. 6 , p. 688.) 

In 0, Section 2, it was shown, that, from an analytical point of view, 
the (spatial) /i-transform of the superdiffusion is given by a certain trans- 
formation of the corresponding semilinear operator. This remains the case 
for the space-time //-transform. 

Lemma 3 (H-transform) 

Let X H be defined by \44\ l- Then 

(a) X H is a(L + -Vj + ^f + ^-,aH; D) -superdiffusion. 

(b) In particular, if X is a time-homogeneous (L, (5, a; D)- superdiffusion, 
and H is of the form 

H(x,t) := e- Xct h(x), (45) 
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where A c is the principal eigenvalue of L + (3, and h is a positive 
solution of (L + (3)h = \ c h, then X H is a(L + ■ V, 0, ahe~ Xct ; D)- 
superdiffusion. 

Remark 6 (Unbounded /3's) 

As it is already the case with the spatial /i-transform for superdiffu- 
sions, it is possible that the coefficient (3 transforms into a new coeffi- 
cient that is no longer bounded. In fact this can be the very definition 
of superdiffusions for certain unbounded 0's (see 0, Section 2 for ex- 
planation). 

Proof of Lemma |3J In order to avoid minor technical inconveniences, 
we implement the method in [g] (see the second paragraph on p. 689.). 
Namely, we use that the Laplace transition functional restricted to the family 
of measures Ai c (D) and the family of functions C^~(D) uniquely determines 
a measure- valued Markov process, and we choose working with these smaller 
spaces rather than replacing Aif(D) and C^(D) by //-dependent spaces. 

Pick v G M C {D), and / G C C {D). Define ^ := v/H(-,s) G M C (D), 
and #(*)(•) := H(-,t)f(-) G C+(D). Obviously, 



exp-(X t ,g^) 



[exp-(X t H ,f)]=W (s) 
By the log-Laplace equation ()40|) . we can continue with 

= exp - ( — - r ,u(; s; t, 

Consider the operator 

A : C 2 ^{D) x C 1 ' ri {R + 

defined by 



(46) 



C V (D) x C^{1 



A(u) := d s u + (L + (3)u — au 2 . 
Define the //-transformed operator A H in the usual way: 



(47) 



A H (u) := —A(Hu). 



H 



Then a direct computation gives 



A H {u) 



VH ~ LH 

—u + o s u + Lu + a—— -Vu + (3u + —^-u 
n ri n 



aHu 



(48) 



(49) 
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Another, trivial computation yields that if 



w( v ;t,/) :=u(;-;t,H(;t)f)/H(;t), 

then v(-, •; t, f) is the solution in (|41l) with .A replaced by A H , and with the 
property in Definition Hfii). Thus the quadruple (L, (3, a; D) transforms into 
the quadruple given in part (a). 

Part (b) is straightforward computation. □ 
Remark 7 

It is precisely equation (|48|) that justifies the name '//-transform'; 
the transformation on the semilinear operator works the same way as 
Doob's ^-transform would work on a linear operator. 
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